NASA/CR— 1999-209313 



Modeling of Flow Transition Using an 
Intermittency Transport Equation 

Y.B. Suzen and P.G. Huang 

University of Kentucky, Lexington, Kentucky 


Prepared under Cooperative Agreement NCC3-590 


National Aeronautics and 
Space Administration 


Glenn Research Center 


September 1999 



Acknowledgments 


This work is supported by NASA Glenn Research Center under grant NCC3-590. The project is part of the Low 
Pressure Turbine Flow Physics program of NASA-Glenn. The contract monitor is Dr. David Ashpis. The authors 
would like to thank Prof. Peter Bradshaw of Stanford University for his review and his comments. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 
Price Code: A03 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 



Modeling of Flow Transition Using an Intermittency Transport 

Equation 


Y. B. Suzen and P. C. Huang 
Department Mechanical Engineering 
Cniversity of Kentucky 
Lexington. Kentucky 10506 


Abstract 


A new trans])ort equation for intermittency factor is proposed to model transitional 
flows. The intermittent behavior of the transitional flows is incorporated into the com- 
putations bv modifying the eddy viscosity, /q, obtainable from a turbulence model, with 
the intermittency factor. //* — y/q. In this paper. Mentor’s SS I model (Menter. 1994) 
is employed to compute/// and other turbulent quantities. 

The proposed intermittency transport equation can be considered as a blending of 
two models Steelant and Dick (1996) and (’ho and Cluing (1992). The former was 
proposed for near-wall flows and was designed to reproduce the streamwise variation of 
the intermittency factor in the t ransition zone following Dhawan and Xarasimha correla- 
tion (Dhawan and Xarasimha, 1958) and the latter was proposed for free shear flows and 
was used to provide a realistic cross-stream variation of the intermittency profile. 

The new model was used to predict the 15 series experiments assembled by Sav- 
ill (1995a, 1995b) including flows with different freest ream turbulence intensities and two 
pressure-gradient cases. For all test cases good agreements between the computed results 
and the experimental data are observed. 
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1 Introduction 


To a large extent, flows in low-pressure turbine applications are transitional, and the 
behavior of the flow transition is strongly affected by the freestrearn turbulence level and 
pressure gradients in the turbine blade passage. In order to predict the losses and heat 
transfer one must be able to predict accurately the boundary-layer development and its 
interactions with flow transition and separation. 

Although a lot of work has been reported on improving the accuracy of the CFD tools 
in engineering predictions, specially in the turbulence modeling areas, CFD prediction 
of transitional flows still remains an almost untouched territory, as compared with other 
areas of CFD developments. A common way to simulate laminar to turbulent transition 
in computations is to switch on the turbulence model (or turbulent eddy viscosity) at an 
experimentally pre-determined transition location. This method is ad hoc and ignores 
the transition physics and the importance of the transitional zone completely. Especially 
for flows where the transitional region covers a large portion of the flowfield, as observed 
in many low-pressure turbine experiments, this practice can lead to severe errors in the 
solution. 

Savill (1993a, 1993b). and West in and Henkes (1997) have tested a large variety of 
turbulence models and compared models' performances in predicting a few of T3-series 
transition flow experiments. They showed that no model could predict both transition 
location and length for a range of flow conditions and concluded that the existing models 
were inadequate to predict flow transition. 

An alternative to this approach is to use the concept of intermittency to blend the flow 
from the laminar to the turbulent regions. This approach, although highly empirical, has 
shown some successes in predicting transition behavior. Dhawan and Narasimha (1958) 
correlated the experimental data and proposed a generalized intermittency distribution 
function across flow transition. The correlation was later improved by Gostelow et 
al. (1994) for flows with pressure gradients subject to a range of freestrearn turbulence 
intensities. 

Solomon et al. (1995), following the work of Chen and T hyson (1971), developed an 
improved method to predict transitional flows involving changes in pressure gradients. In 
this model, the effects of changing streamwise pressure gradient on the breakdown physics 
and spot spreading rates are taken into account. This is accomplished by varying the spot 
spreading angle and propagation parameter through the transition zone according to the 
local pressure gradient parameter. 

Steelant and Dick (1996) proposed a transport equation for intermittency, in which 
the source term of the equation is developed such that the 7 distribution of Dhawan 
and Narasimha (1958) across the transition region can be reproduced. Steelant and 
Dick used their model, coupled with two sets of conditioned Navier-Stokes equations, to 
predict transitional flows with zero, favorable, and adverse pressure gradients. However, 
since their technique involved the solution of two sets of strongly coupled equations, the 
method is not compatible with existing CFD codes, in which only one set of Navier-Stokes 
equations is involved. Moreover, the model was designed to provide a realistic streamwise 
7 behavior but with 110 consideration of the variation of 7 in the cross-stream direction. 
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Clio and ('hung ( 1992) developed a A*— e — y turbulence model for free shear flows. Their 
turbulence model explicitly incorporates the intermit tency effect into the conventional A “— ( 
model equations by introducing an additional transport equation for "} . They applied this 
model to compute a plane jet, a round jet, a plane far wake, and a plane mixing layer with 
good agreements. Although this method was not designed to reproduce flow transition it 
provided a realistic profile of y in the cross-stream direction. 

In the current paper, a new t ransport equal ion for intermittency factor is proposed. 
The main motivation in the development of the new model is to combine the best features 
of the existing transition models in the proposed model. The model can not only reproduce 
the intermittency distribution of Dhawan and Narasimha in the streamwise direction but 
also is able to provide a proper variation of y in the cross-stream direction. Our aim is 
to propose a model that will predict flow transition under the influences of freest ream 
turbulence and pressure gradients. 

A review of a number of transition models used in the current study is given in sec- 
tion 2. In section 3, the new model, combining the best features of the models discussed 
in section 2. is proposed. Section 1 discusses the numerical issues regarding the imple- 
mentation of the current method in the existing CFI) codes, including our choice of the 
correlations to predict the onset of the transition under the influences of freest ream tur- 
bulence and pressure gradients. Section 5 shows the comparisons of the new transition 
model against T3 series experiments of Savill (1993a, 1993b). Concluding remarks are 
provided in section G. 


2 Transition Models 

2.1 Dhawan and Narasimha 

Dhawan and Narasimha (1958) correlated a range of experimental data and provided the 
following expression for streamwise intermittency distribution; 


I - exp[-(.r - .r t ) 2 ncr/r] (.v > .r,) 
0 U<.r,) 


where ( is the freest ream velocity; .r, is the point of transition onset: n is the spot 
formation rate and rr is the spot propagation parameter. 

There are two ways to express equation (I) in dimensionless form: one way. as given 
by Narasimha (1985), is to express equation (1 ) in terms of a nondimensional breakdown 
parameter, V = nrrQf/is, where 0 t is the momentum-thickness at the onset point of 
transition, .r — .r f , and v is the kinematic viscosity: and the other way is in terms of 
dimensionless spot formation rate, u ~ nv 2 /C r3 , as used by Mayle (1991). Hence, the first 
part of equation (1) can be written either as: 


or 


5 = 1- exp 



2 A' 


/ = 1 ” exp 



2 

naRel 


( 3 ) 
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Gostelow et al. (1994) conducted measurements of boundary layer transition for six 
different turbulence levels and a wide range of pressure gradients and correlated the 
nondimensional breakdown parameter A as; 


V = 


0.86 x 10~ 3 exp[2.134Afl \n{Tu) - 59.23A* - 0.564 ln(7u)] \ 0 < 0 
0.86 x 10 -3 exp[— 0.564 ln(T»)] X 9 > 0 


(4) 


where the pressure gradient parameter. \g = ( 9 2 /v)(dU/dx ), and the freestream turbu- 
lence intensity, Tu, are evaluated at the onset point of transition. 

In contrast, for zero-pressure-gradient flows, Mayle (1991) correlated the spot produc- 
tion rate using intermittency measurements as; 


her = 1.25 x 10 7V 


r/4 


(5) 


The constant in equation (5) may vary slightly depending on the choice of the data: Mayle 
also reported a value of 1.5 x 10~ n in the same paper. When subject to pressure gradients. 
Mayle (1991) normalized her by its value at zero-pressure-gradient, ( ha)zpa , and showed 
a plot of her /(hcr)zpG versus the acceleration parameter, A = (u/U' 2 )(dU/dx) over a 
wide range of freestream intensity, Tu. Steelant and Dick (1996) proposed the following 
correlation to fit Mavle’s data; 


ner 


{ner )zpg 

A comparison of equations (2) am 


_ J (4747hr 2 - 9 ) 1-exp(2xl ° 6/v) , A < 0 

3227A' 0 59sr ' 


10 - 


I< > 0 


( 6 ) 


[3) indicates that: 
A = herR(g t 


If both Gostelow et al. (1994) and Mayle (1991) used the same data to correlate their 
dimensionless quantities. A and her, respectively, an equation for the Reynolds number at 
the onset of transition can be identified. For example, for a zero- pressure-gradient case. 
Gostelow et al.'s (1994) correlation, equation (4). gives; 


A 


-0.564 


0.86 X 10" 3 7V 
5) into equation (7). yields; 


Substituting equations (8) and 

Re 0 , = 410 r<r°- 771 

Since equation (9) is different from Mayle’s correlation (Mayle, 1991); 

Rio, = 4007’« -u ' 625 


( 8 ) 


(9) 


( 10 ) 


nor does it resemble the well-known Abu-Ghannam and Shaw correlation for zero-pressure- 
gradient flows (Abu-Ghannam and Shaw, 1980); 


Rig, - 163 +exp(6.91 — Tu) 


11 


the two correlations, equations (5) and (8), are not identical. We concluded that the 
choice of the transition correlations using either N or her must be tested in company with 
a proper choice of the correlations for the onset of transition. 
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2.2 Solomon et al. 


Typic ally for flows in low-pressure turbines, the onset of transition and the end of transi- 
tion are in different pressure gradient regions, and to derive transition correlations using 
only parameters at the onset point of transition may result in an incorrect prediction 
of the transition length. Following the work of Chen and Thyson (J971). Solomon et 
al. (1995) developed a modified method for calculating intermittencv in transitional flows 
with changing pressure gradients. 

The model of Solomon et al. (1995) accounts for the effects of changing streamwise 
pressure gradient on the breakdown physics and spreading rate of turbulent spots by allow- 
ing the model constants to vary continuously with the local pressure gradient parameter, 
A*. While the spot generation rate is assumed to depend only on the local conditions at 
the transition onset point using the relation of Costelow et al. (1994), equation (4). the 
spot spreading half-angle o and the spot propagation parameter cr are given by: 

^ = 0.03 + (0.37/(0.48 + 3.0exp(52.9A*))) (1*2) 


o = 4 + (22.14/(0.79 + 2.72exp(47.63A*))) 


(13) 


The intermittencv distribution, retaining the concentrated breakdown hypothesis of 
Narasimha. is given by: 


y = 1 — exp 


r J i a d.v r 

r t 

n / r . / 

tan odv 

Jx t tano ( Ja 

t 


(14) 


where* the spot generation rate. n. is calculated using A from equation (4). 

Solomon et al. (1995) tested their model for several pressure gradient flows and showed 
improved agreement with the experimental data. 


2.3 Steelant and Dick 


Steelant and Dick (1996) developed a transport model for intermittencv to be used in 
conjunction with conditioned Navier-Stokes equations. They derived a model equation 
starting from the intermit tency distribution of Dhawan and Narasimha (1958). Differ- 
entiating equation (1) along the streamline direction, .s. they arrived at the following 
transport equation: 


Of 


Vpin 

+ ().V 


+ 


dpn 

(h) 


(1—5 )p\/ ii 2 T e 2 T(s) 


(15) 


with 


T(.s) = 2/(.s)/'(.s) 


(16) 


The function d(.s) represents the (ncr/C)(.r — .r f ) term in equation (1). flic function f(s) 
is formulated to account for distributed-breakdown and is given as: 


/(*) 


».S M + 6.s /:i + cs' 2 + <i + e 

gs rA -)- h 


(17) 
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where the coefficients are; 


a = b = -0.4906 c = 0.204(^)" 0 - 5 


d = 0.0 e = 0.04444(^)- ]r ’ g = 1.0 It = 10 c 


Equations (5) and ( 6 ) were recommended to evaluate the values of ncr. In the above equa- 
tions* l is the freest ream velocity at the transition location and the streamline coordinate. 
<$. is defined as; 


.s 


udx + vdy 

\fu 2 + V 2 


( 18 ) 


and = $ — $ f , where $ f is the transition location. 

Steelant and Dick (1996) tested their model in conjunction with two proposed sets 
of conditioned averaged Navier-Stokes equations for zero, adverse, and favorable pressure 
gradient flows and their results showed that the model performed well for all cases. While 
the intermittencv transport equation reproduces the streamwise intermittency distribu- 
tion of Dhawan and Narasimha, it gives rise to a fairly uniform y distribution in the 
cross-stream direction. This is somewhat inconsistent with the experimental observation 
of Klebanoff (1955), in which the variation of intermittency in the cross-stream direction 
has been proposed empirically as; 


7(y) = 1[1 - erf(C)] (19) 

with 

C = ' r > [ l(|: ) - 0 ‘S] (20) 

where S" is the displacement thickness. 


2.4 Clio and Chung 

Cho and Chung (1992) developed a h — t — 7 turbulence model for free shear flows. In their 
model the intermittency effect is incorporated into the conventional k — t turbulence model 
equations explicitly by introducing a transport equation for the intermittency factor 7 . 
The eddy viscosity is expressed in terms of k, c. and 7 . Their intermittency equation was 
given as; 


'dxj 


D, + 5, 


( 21 : 


The diffusion term, D~,, is represented by; 


£>, = 


d 


(1-7) 


i't d'j 


d.Vj L <7y d.Vj . 


•») 
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and the source term. , is; 


V = r 9 ,7(i -7) 


f\, + I\,n , „ k 2 ch <h 


k 


+ ( 'q>' 


( dxj d.Vj 


C.37(i - 7) r r 


( 23 ) 


vv 


here 


l\-.s = ~ u i u jTT~ ^ j) 

dxj 


(24) 


dui 

‘k.n = Hi'U/Ti (*=./) 

■ dr; 


(25) 


and 


r = 


k :, k 2 u, da, (h t 
( 2 (»a»a) 1/2 dxj i)Xj 


(26) 


The modeling constants are; 


rr S) =1.0 ( ' a i = 1.6 ( >2 =0-15 ( V- ( = 0.16 

Clio and Chung (1!M)2) tested t heir model for a plane jet., a round jet. a plane far wake, and 
a plane mixing layer and showed improved model performance. Although the model dot's 
indeed produce a very realistic y profile for turbulent free shear Hows, it is not designed 
to predict flow transition to t urbulence. 


3 New Transport Model for the Intermittency 

A new transport equation for intermittency is proposed. The main objective of the new 
model is to reproduce the intermittency distribution of Dliawan and Narasirnha in the 
st ream wise direction and at the same time to give rise to a realistic variation of the 
intermittency in the cross-stream direction. In order to accomplish this, a transport 
model for intermittency, blending Steelant and Dick's and (To and Chung's models, is 
proposed. 

The production term for the new model is a mix of the generation terms of Steelant 
and Dick and of Cho and ( Tiling. The first term, 7C is from St cedant and Dick, aiming 
to reproduce the intermittency distribution of Dliawan and Narasirnha. The formula for 
T ( , is given bv: 

To = Cops/ a 2 + rO(.s) (27) 


where 


■*(") = 2/(.s)/'(.s) 
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One major aspect in which term T 0 differs from Steelant and Dick’s model is the distributed- 
breakdown function f(s). Steelant and Dick’s model was calibrated for use with condi- 
tioned Navier-Stokes equations and when used with the current approach it allows ex- 
cessively long regions of distributed-breakdown. In order to a obtain shorter distance for 
distributed-breakdown the function f{$) is modified in the new model. 

The function f(s) used in the new model has the same form as Steelant and Dick's 
model; 


a*' 4 + bs a + c.s /2 + ds' + e 

/(s) = ^TT, 

where the coefficients are; 

ill (7 , 77(7 a c; 

a = 50^— 6 = -0.4906 c = 0.204( — ) 

d = 0.0 c = 0.04444 ( jr)" 1 ' 5 h = 10c g = 50 


(29) 


(30) 


These coefficients are the same as those used in Steelant and Dick's model except for a and 
g. As will be seen in the comparisons of the models shown in section 5, the use of equation 
(17) together with the current intermittency approach, which makes use of intermittency 
to modify turbulent viscosity directly (see section 4 for details), has a tendency to delay 
the onset of flow transition. Since the current approach does not use the conditioned 
Navier-Stokes equations, this adjustment of coefficients, leading to a shorter distance 
for the distributed-breakdown function f{s), enables a faster response of the mean flow 
solutions to the intermittent effects. 

Two major production terms from Cho and Chung's model are T\ and 75- These two 
terms are used in the form {T\ — 75) in the model. The term T] mimics the production 
of turbulent kinetic energy, P/., and is given by; 


T-r ^_C l7 dm 

Tl ~ Cll T ~ i r^dTj 


with the shear stresses defined as; 


(31) 


T-,., = (h 


dm dxij 2 dtik 

dx, dx, ndxk' JlJ 


6 ,, 




(32) 


The term T 2 represents the production resulting from the interaction between the mean 
velocity and the intermittency field and is given by; 


T 2 = C 2 HP- 




u, 


duj d' s 


(«*«*) 1/2 dx j dx, 


(33) 


The production terms 7’ 0 and ( 7j — T > ) are blended by using a function F to facilitate 
a gradual switching from Steelant and Dick's To to Cho and Chung's (Tj — T 2 ) inside the 
transition region; 


l\ = ( 1 — F)T q + F(7j — T 2 ) 


(34) 
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A nondimensional parameter. k/Su, is chosen to construct the blending function f\ 
where A* is the turbulent kinetic energy and S is the magnitude of the strain rate. This 
parameter increases rapidly with distance away from the wall inside the transition region. 
Using the solution obtained bv Dliawan and Narasimha's model coupled with artificial 
cross-stream profiles provided by KlebanofTs formula, equation ( 19 ). one may divide the 
transition zone into two regions by taking a diagonal cut between the point at the edge 
of the boundary layer located at the beginning of transition (7 = 0) to the point at the 
wall located at the end of transition region (7 = 1). T he relation between 7 and {kjSu) 
along this diagonal line is then approximated by t he following correlation; 

k/Su — 200( 1 — 7 01 )°' 3 ( 35 ) 


This line corresponds to the border between the Steelant and Dick and the (T10 and 
(Tiling models: i.e. below this line, 7 ’o is active and above this line ( / 1 — 1 >) is active. 
In order to facilitate a gradual switching from 7 0 to (T\ — 7 V). the following blending 
function is proposed: 


F 


t anh 4 


k I Su 

200 ( 1 - 7 0 - 1 


( 36 ) 


As can be seen from equation ( 36 ). when k/Su >> 200(1 — ^ u l )°* 3 (a position above the 
cutoff line). F = 1 and the model switches to (T10 and Chungs model; when k/Su << 
200(1 — 7°‘ 1 ) u '’ 3 (a position below the cutoff line), F — 0 and the model becomes Steelant 
and Dick's model. Outside the transition zone, the model switch to the Clio and (Tiling s 
model except for the very thin region close the wall. 

In figure 1. the [imposed correlation for the cutoff line, equation ( 35 ), is compared 
with data extracted along the cutoff line, defined by drawing a straight line between the 
point at the edge of the boundary layer located at the beginning of transition to the point 
at the wall located at the end of transition region, of the T 3 A solution obtained using the 
[imposed transition model. As can be seen from the figure, the definition of the cutoff 
line [imposed in equation ( 35 ) is satisfactory. 

An additional diffusion-related production term is introduced by Cho and (Tiling as: 


7 T - C V 


k 2 c) 7 rh 

( c)xj c)x J 


( 37 ) 


This term is kept active over the entire flowfield. that is. 110 blending is applied to this 
term. 

Diffusion of 7 is represented by the following term: 


= jjk j[(l - 7 ) 7 ( 7 .,,// + ( 1 - 7 K,/'/] 


-I 

J 


The final form of the model is: 

0/ n , Optij'i 


rjl + O.r j 


= (1 


(I - F)T 0 +F( 7 \-T 2 ) 


+ 7 :i + D- 


(:W) 
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or 


dpj_ dpujj _ 

dt dx, 1 


(1 - F)C 0 pVu 2 + r' 2 T(.s) 


, r /Cn dm 

+ F (— ^57" 
^ r ** a ‘ ^ 

"f (-3 P 


- C 2 lp 


k 3 /' 2 Uj dll, d~ 
t {u k u. k ) x / 2 c)xj dxj j 


f dxjdxj 

+ ^-((( 1 - 7)7 v^p + ( 1 - 7 ) (7 ~>tPt)^r) 

where the modeling constants are. 

<t-v ( = oy, = 1.0 Co = 1.0 t'i = 1.6 C -2 — 0.16 ('?, — 0.15 


( 40 ) 


4 Implementation of the Transition Model 

The intermit tency concept can be incorporated into the computations either by using 
conditioned-average Navier-Stokes equations (Steelant and Dick, 1996: Libby. 1975) or 
simply by multiplying the eddy viscosity obtained from a turbulence model, /q, bv the 
intermittency factor. 7 (Simon and Stephens, 1991). The major difficulty of the former 
method is the requirement to solve two sets of highly coupled conditioned Navier-Stokes 
equations. This method is computationally expensive and is not compatible with existing 
CFD approaches. The latter approach is recommended in the current paper. Simon 
and Stephens (1991) showed that by combining the two sets of conditioned Navier-Stokes 
equations and making the assumption that the Reynolds stresses in the nonturbulent part 
are negligible, the intermittency can be incorporated into the computations by using the 
eddy viscosity, //*. which is obtained by multiplying the eddy viscosity from a turbulence 
model, ftt, with the intermittency factor. 7. That is, //* = 7 /q is used in the mean flow 
equations. 

To allow the intermittency to have full control of the transitional behavior, the turbu- 
lence model selected to obtain /q must produce fully turbulent feature before transition 
location. Menter's SSL model (Menter, 1994) can he shown to produce fully turbulent 
flow in the leading edge of the boundary layer. It is therefore recommended to be used 
as a baseline model to compute /q and other turbulent quantities in the calculations. 
Menter's SST model is provided in the Appendix. 

The value of no used in evaluating the constants given by (30) is provided by the 
M a vie correlation; 

ho = 1.8 x 10“ 11 Tu‘/ 4 (41) 

It should be noted that a value of 1.8 x 10 -11 was used in the current work to give a slightly 
better lit of Mayle's data. When flows are subject to pressure gradients, the correlation 
of Steelant and Dick, equat ion (6), is used. 

The current intermittency approach was applied in conjunction with the correlation 
of Huang and Xiong (1998) for the onset of transition: 

Re et = ( 120 + 150r»- 2/3 )coth[4(0.:f - I\, x 10 r, )j (42) 
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where l\ f is the minimum value of the acceleration parameter in the downstream decel- 
eration region. This correlation was found to provide a slightly better approximation 
than the correlation of Abu-Gliannam and Shaw (1980) in favorable pressure gradients, 
as discussed by Huang and Xiong (1998). 

The computations are performed with a boundary layer code which solves the mean 
flow, turbulence model and intermit tency equations using second-order finite-volume 
method. In the computations. 175 grid points, expanding from the wall to the freest. ream, 
were used in the cross-stream direction for all cases. The values for the first point away 
from the wall were kept between 0.1 and 0.15 for all cases. The solutions were obtained 
by using 1000 streamwise steps for all cases. This corresponds to maximum dimensionless 
streamwise step sizes, A.r + . of 49. 68, 47, and 47 for cases I 4A. 14 B, T4(T, and T4C2 
respectively. These step sizes and cross-stream grid points were found satisfactory by 
performing a careful grid-independency check, in which the step sizes and grid spacing 
were both decreased by half: no effect on the solutions was found. 

At the inflow, a top-hat velocity profile is prescribed. Inlet turbulent kinetic energy 
is fixed according to the freestream turbulence levels and the energy dissipation rate is 
adjusted according to the decay of the freestream turbulence, as will be illustrated in the 
next section. 


5 Results and Discussion 

The new transition model is used to predict the experimental test cases assembled by 
Sax ill (1994a, 1994b): T4A, T4U, T4C1, and I 4(42. These experiments were specially 
designed to test the ability of turbulence' models to predict the effects of freestream tur- 
bulence on the development and subsequent, transition of a laminar boundary layer under 
zero and varying pressure gradient condit ions. Gases T4A. I 4B are zero-pressure-gradient 
flows and 4 4(41. 14(42 are cases with continuous change in pressure gradient represent- 
ing an aft-loaded turbine blade. Comparisons were performed for these cases with the 
new transition model against the models of Dhawati and Narasimha. Solomon et al.. and 
Steelant and Dick. 

In all computations, the inlet conditions were calibrated to match the experimental 
decay of turbulence using the SST model. In each case, the inlet turbulent kinetic en- 
ergy was fixed by the experimental freestream turbulence level and the matching of the 
freestream turbulence decay provided an estimated value for the dissipation rate of tur- 
bulent kinetic energy. (. (or the value of /q///) at the inlet. Once the value of < (or //*///) 
was determined this value was used in subsequent transition computations. 

In all cases, the onset of transition was specified according to the correlation of Huang 
and Xiong (1998). equation (42). The data of Savill were found to correlate well with equa- 
tion (42) considering the uncertainties in experimental data; see Huang and Xiong ( 1998). 
For 4 4 A, T4B. T4(\l. and T4C2. the correlation gave the transition locations as Ii($ t = 
257, 204. 194. and 297, respectively, while the experimental data showed Rts t — 272. 182. 
211, and 478, respectively. The determination of the experimental onset of transition is 
rather arbitrary. In our mode] calculations, we deliberately used the correlation to predict 
the onset of transition in order to reduce the dependency on experimental inputs. 
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T3A experiment of Savill (1993a) corresponds to a zero-pressure-gradient flow over a 
flat plate at Re = 3.6 x 10 5 per meter. The freestream turbulence intensity at the leading 
edge of the flat plate is 3.35%. The decay of freestream turbulence intensity is matched 
with the experimental data by specifying /<,//< = 7.6 at the inlet, as shown in Figure 2. 

The computed skin-friction coefficient distributions are compared with the experi- 
mental data in Figure 3. Since T3A case is a zero-pressure-gradient case, Dhawan and 
Narasimha’s and Solomon et al.'s models are identical. Steelant and Dick’s model showed 
that the transition was somewhat delayed and the length of transition was slightly wider 
than the experimental data when compared with the data and the other models. In 
contrast, the current model predicted the length of the transition region well and also 
displayed an overshoot of skin-friction coefficient, Cj. at the end of transition region. 
It should be noted that the failure of Steelant and Dick's model is because the original 
distributed-breakdown function, equation (17). was calibrated against the conditioned 
Xavier- Stokes method. To allow' for a. faster repsonse to flow transition when coupled 
with the current approach, a modified distributed-breakdown function, equation (29). is 
recommended . 

The comparison of the Reynolds number based on momentum-thickness, Reg. for T3A 
case is shown in Figure 4. As can be seen from this figure, the result of the new model is 
also in good agreement with the experimental data. 

One of the major features of the current model is its ability to reproduce realistic cross- 
stream intermit tency profiles. The predicted intermit tency profiles at various stream wise 
stations through the transition zone are shown in Figure 5. The profiles exhibit a peak 
between y/S* = 1 and y/S* = 2 then drop off toward zero near the edge of the bound- 
ary layer, around y/S* = 8. These features are consistent with the trends observed in 
experimental data of Solin and Reshot ko (1991), and Gostelow and Walker (1991). 

The current intermittency model does not affect the ability of the turbulence models 
in the fully developed region, as can be demonstrated in Figure 6. In figures 6(a) and 
6(b). the velocity profiles obtained using the current intermittency model are compared 
with those predicted using the SST model without the intermittency modification for 
Re $ = 1,000 and 5,000. respectively. For Reg = 1000. the velocity profiles differ slightly 
near the freestream, indicating the Cj value predicted bv the current model is slightly 
larger. This larger value of the skin friction is associated with the overshoot of the skin 
friction near the end of the transition region, as shown in Figure 3. Further downstream, 
at Reg = 5000 the velocity profile obtained from the new model is essentially coincident 
with the one obtained by using the SST model alone, as shown in figure 6(b). These results 
clearly show' that the new intermittency model does not affect “'good solution behavior 
of the t urbulence model the in the fully developed turbulent region. 

The second test case is the T3B case of Savill (1993a). T3B case is also a zero- 
pressure-gradient flow with a freestream turbulence intensity of 6% at the leading edge of 
the flat plate. The freestream Reynolds number is 6.3 x 10 r> per meter. In order to match 
the experimental decay of free stream turbulence intensity, as shown in Figure 7, /q//t is 
specified as 60 at the inlet. 

The skin-friction coefficient distributions obtained by the models are compared with 
the experimental data, as shown in Figure 8. The transition length predicted by the 
Dhawan and Narasiinha correlation is longer than the data and those predicted by the 


NAS A/CR— 1 999-2093 1 3 


12 



other two models. Steelant and Dick’s model showed a delay of transition because of 
the use of the distributed-breakdown function, equation (17). The new model produces 
the shortest transition length, and is closest to the experimental data. The discrepancy 
between the current solution and the experimental data for this case is partly due to the 
specified transition onset point. 

In Figure 9 . the momentum-thickness Reynolds number is compared to the experimen- 
tal data. The solution of the new model is the closest to the experimental data followed 
bv Steelant and Dick's model. 

The 7 profiles from the new model at various streamwise stations through the transi- 
tion zone are shown in Figure 10. The profiles have the same characteristics as the T3A 
case, but due to high free stream turbulence intensity, the results showed that the peaks in 
the profiles are less pronounced for this case and the spread of the intermit tency appears 
to be wider across the transition region. 

The next cases are the T3(T and T3C2 experiments of Savill ( 1993a). These cases are 
variable pressure gradient flows representing actual turbine characteristics. Both cases 
have similar pressure variations, but with different Reynolds numbers and fret' stream 
turbulence intensities. The pressure coefficient distributions along the Hat plate for these 
cases are shown in Figure 11 . In the figure, the symbols are the experimental data points 
and the lines are the curve fitting to the data points. These ( ' () profiles were input in the 
boundary layer code to perform the calculations. 

T3C1 case has an inlet Reynolds number of 4.1 x 10 5 per meter. Turbulence intensity 
at the leading edge of the plate is 7.78 ( /( and a value of fit /ft — 30 is used at the inlet in 
order to match the decay of free stream turbulence intensity, as shown in Figure 12. 

A comparison of skin-friction is shown in Figure 13 for T3CT case. All models except 
the current, model overpredicted the transition length. The solutions of Dhawan and 
Narasimha model and of the model of Solomon et al. are slightly different with the latter 
predicting a longer transition length. Again. Steelant and Dick’s model showed a delay of 
the transition due to a slower distributed-breakdown function, equation (17), being used. 
Overall, the new model showed a fairly good comparison of the transition behavior. 

The lit o distribution along the flat plate is compared to experimental data in Figure 11 
for the T3(T case. In general, the new model reproduced the development of the boundary 
layer better than all other models. 

The intermitt ency factor profiles from the new model at various streamwise stations 
through the transition zone are shown in Figure 15 for the T3(M case. The characteristics 
of the profiles are similar to the T3B cases with mild peaks showing around y/S m — 1.0 
and the profiles decaying to zero around y/$ M — 10 , near the edge of the boundary layer. 

T3C2 case has an inlet Reynolds number of 3.5 x 1 G :1 per meter and a free stream 
turbulence intensity of 2.8 ( X at the leading. The decay of free stream turbulence intensity 
was matched with the the experimental data by specifying ////// = 5 at the inlet, as shown 
in Figure 10. 

The predicted skin-friction coefficient distributions are compared with experimental 
data in Figure 17. Both the models of Dhawan and Narasimha and of Solomon et al. 
give rise to a too early transition behavior and also a longer transition length. On the 
other hand. Steelant and Dick’s model did not predict the onset of transition upto the 
end of computational domain. Again this is caused by the original distributed-breakdown 
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function being used, as discussed earlier. Overall, the current model shows a fairly good 
comparison with the data. 

The Reynolds number based on momentum-thickness along the flat plate is compared 
to experimental data in Figure 18. Overall, the new model reproduces the development 
of boundary layer fairly well. 

The intermittency factor 7 profiles at various stations through the transition zone are 
shown in Figure 19 for the T3C2 case. It should be noted that T3C2 has a lower freestream 
turbulence intensity than T3C1. As a results, the profiles show more pronounced peaks 
between y/S* = 1 and y/8* = 2 and the spread of the intermittency is less wider across the 
transition zone. These features are similar to the comparison of the zero-pressure-gradient 
cases, T3A and T3B. 


6 Concluding Remarks 

A new transport equation for intermittency is developed for modeling transitional flows 
including influences of free stream turbulence and pressure gradients. The model can 
be considered as a mix of two existing models, aiming to reproduce the intermittency 
distribution of Dhawan and Narasimha (1958) in the streamwise direction and to produce 
a realistic variation of intermittency in the cross-stream direction. 

The new model was tested for both zero-pressure-gradient and variable pressure gra- 
dient flows with different freestream turbulence intensities. The new model was used in 
conjunction with Mentor's SST model to predict the T3A, T3B, T3C1. and T 3 C 2 ex- 
periments of S a v ill (1993a. 1993b). Comparisons of Cj distributions, the development of 
the boundary layer, and intermittency profiles were made for all cases. The new model 
showed a good transition behavior for all cases. 
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Appendix — Menter’s SST Model 

Menter's SST model (Menter, 1994) is based on a mix of two-equation k — *; and k — < 
turlndence models using a blending function F\. The model can be written as: 
A-eqnation: 


dpk Opujk 


dt 




d dk 

Pi. - 0.09 puik + —((/' + Vk/h) — 

d.Vj d.r , 


(A-J 


-equation: 


dpo: doii.u; c ,, 

+ - L 7 — = ~Pk- 


dt 


d.r 


ip 


, t) d^' 

dp^‘ + — ((// + — 

d.r , d.r , 


+ 2/r( 1 — F\ )o jJ 2- 


1 dk cL; 

d.r j d.r j 


(A-2) 


The constants, c. A, cr k and a^. are given by the following general expression: 


o — I' \ o i + ( 1 — F\ )&2 ( A •! ) 

where o represents any one of these constants; Oi represents any constant in the k — u.’ 
model; o 2 represents the corresponding constant in the k — ( model. These constants arc 
defined as: 

Set 1 ( k — u,’ ) : 

cr, ;] = 0.85 ^, = 0.0 A, = 0.075 r, = 0.552 (A-4) 

Set 2 (k — ( ): 


a ) i: 2 = 1 .0 = 0.856 d> = 0.0828 c 2 = 0.44 

The production term is given as: 


Pk = 


T u 


dll; 

d*j 


where 


(A-5) 


(A-6) 


Tij = dt 


duj 

a 7 - + 


ditj 

d.r t 


2 du k 

3 d.r, 



The blending function l\ is defined as: 


(A-7) 


P\ = tanh(arg‘| ) 


(A -8) 
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with 


argj = min 


max 


500// 
0.09u></ d 2 u.’ 


4/?^2 1 


(A-9) 


where <! is the distance to the closest wall and f is the positive portion of the cross- 
diffusion term in Equation (A-2). 


CDkuj = max 


lP<7* 2 - 


$A- du: 
dxj Ox, 


: 10 


-20 


( A - 1 0 ) 


The kinematic eddy viscosity is defined as; 


i't = 


a\k 

max(fl|u.’; 0 F 2 ) 


(All) 


where ft is the magnitude of vorticity and a x = 0.31. The function F 2 is given hy: 


with 


F 2 = tanh(arg^) 


arg 2 


2 y/k 500 v 
0.09u.y/' d 2 ^ 


(A- 12) 


(A-13) 
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Equation (34) 
Computational Result 









Figure o: Iiit.erniittencv factor profiles for T3A cas<' 
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(a) Comparison of velocity profiles at Re e =1000. 



(b) Comparison of velocity profiles at Re e =5000. 


Figure 6: Comparison of velocity profiles at two stations. 
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Figure 7: Comparison of freest- ream t urbulence intensity for T 3B case, 



Figure 8: Comparison of skin-friction coefficient for T3B case. 
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